Non-Adiabatic Holonomic Quantum Computation 
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We develop a non-adiabatic generalization of holonomic quantum computation in which high- 
speed universal quantum gates can be realized by using non-Abelian geometric phases. We show 
how a set of non-adiabatic holonomic one- and two-qubit gates can be implemented by utilizing 
optical transitions in a generic three-level A configuration. Our scheme opens up for universal 
holonomic quantum computation on qubits characterized by short coherence times. 
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Circuit-based quantum computation relies on the abil- 
ity to perform a universal set of quantum gate opera- 
tions on a set of quantum-mechanical bits (qubits). A 
key challenge to achieve this goal is to find implementa- 
tions of gates that are resilient to certain kinds of errors. 
Holonomic quantum computation (HQC) y| is a general 
procedure to build universal sets of robust gates, by using 
non-Abelian geometric phases. 

Adiabatic holonomic gates have been proposed for 
trapped ions [2], superconducting nanocircuits [3|, and 
semiconductor quantum dots [4]. These gates have 
turned out to be difficult to realize experimentally. An 
obstacle is the long run-time required for the desired 
parametric control associated with adiabatic evolution. 
In other words, as these gates operate slowly compared 
to the dynamical time scale, they become vulnerable to 
open system effects that may lead to loss of coherence. 
On the other hand, if the run-time is decreased in or- 
der to shorten the exposure, non-adiabatic corrections 
start to become significant and the parametric control is 
lost. These problems can be tackled by using Abelian 
non-adiabatic geometric phases [5J to realize quantum 
gates [6]. However, these gates are limited to commuting 
operations and thus cannot perform universal quantum 
computation. 

To combine speed and universality, we propose in 
this Letter a non-adiabatic generalization of HQC. We 
demonstrate an experimentally feasible optical scheme 
to implement a universal set of holonomic one- and two- 
qubit gates for non-adiabatic optical transitions in three- 
level A configurations. 

We first show how quantum holonomy arises in non- 
adiabatic evolution. Consider a quantum system charac- 
terized by an iV-dimensional Hilbert space. Let C : t i-> 
M{t) be a smooth path of X-dimensional subspaces {K < 
N) generated by the Hamiltonian H{t). Mathematically, 
C is a path in the Grassmann manifold G{N; K), being 
the set of complex ii'-planes in C^. If P{t) is a projec- 
tor onto the subspace M{t), then ihP{t) = [H{t),P{t)]. 
We consider t G [0,t] such that P{t) = P(0), i.e., loops 
in G{N;K). Now, let |Cfe(t)), k = \,...,K, be a once 



diffcrcntiable set of orthonormal ordered bases of M{t) 
along C such that |Cfe(''')) = ICfe(O)). The set of all bases 
forms a Stiefel manifold S{N;K), which is a fiber bun- 
dle with G{N; K) as base manifold and with the set of 
K X K unitary matrices as fibers [7| . A lift of the loop 
C in G{N; K) to a loop C in S{N; K) corresponds to a 
single- valued choice of gauge. A gauge transformation is 
a unitary change of bases over C. The final time evolu- 
tion operator projected onto the initial subspace may be 
written as, {h — 1 from now on) t8| 

K 

U{r,0)= J2 (Te^fo^Mt)^"it))dt\ |^,(o))(0(0)|, (1) 
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kl 



where T is time ordering. Here, Af^i{t) = i(Cfe(i)IC/(^)) 
and Hki{t) = {Ck{t)\H{t)\Ciit)) are Hermitian K x K 
matrices. A transforms as A — >■ V' AV + iV^V under 
the gauge transformation \Ck{t)) -> J2i=i\Ci{t))Vik{t), 
where V{t) is a once differentiable family of unitary 
K X K matrices. Thus, A transforms as a proper vector 
potential. The unitary 



U 



(2) 



is the holonomy matrix in the representation associated 
with the loop C, where A-ki — i{Ck{t)\d(i{t)) is the cor- 
responding matrix-valued connection one-form. Note 
that P is path ordering in S{N; K) and that U — > 
V^(0)C/V'(0) under a gauge transformation. This gauge 
covariance essentially means that the holonomy matrix is 
a property of the loop C in the base manifold G{N; K), 
i.e., U=U{C). 

Non-adiabatic HQC is realized under the following two 
conditions: (i) there should exist physically accessible 
loops C : [0, r] 9 t !->• M(i) of subspaces along which the 
Hamiltonian matrix Hi^iit) — {C,k{t)\Pl {t)\C,i{t)) vanishes 
for all pairs |C/c(0)i I0(^)) € M{t); (ii) there should exist 
two such loops C and C", based on the same point in 
G{N;K), for which the corresponding U{C) and U{G') 
do not commute. While the first condition assures that 
the evolution is purely geometric, the second one is nec- 
essary to realize universality. Under conditions (i) and 



(ii), there is a set of quantum gates 



K 

U{r,0) = U{C) = Y, C/,,(C)|a(0))(0(0)| (3) 

kd=l 

that may be able to perform non-adiabaticahy any com- 
putation on qubits encoded in M(0) based purely on the 
geometric properties of G{N;K). We demonstrate that 
these conditions can be met in a generic three-level A 
configuration, by means of which a universal set of one- 
and two-qubit gates can be realized. 
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FIG. 1: (Color online) Setup for non-adiabatic holonomic 
one-qubit gate in a A configuration. A pair of laser pulses 
couple two ground state levels and 1 to an excited state e 
(left panel) . The two ground state levels define a single qubit 
and the laser parameters satisfy \ujo\'^ + \i^if = 1. The dark 
state |d) = —uj-i\Q) -|-a;o|l) is decoupled from the bright state 
|6) = ljq\Q) -\- uil\l) by choosing time-independent cjq and cji 
over the duration of the pulse pair. The system thereby per- 
forms Rabi oscillations between the bright and excited states 
with frequency Q.{t) (right panel). The evolution of the qubit 
subspace is purely geometric and becomes cyclic after com- 
pleting a Rabi oscillation by choosing f2(f) to be a real-valued 
■K pulse. The resulting unitary quantum gate operation act- 
ing on the qubit is determined by the holonomy of the cor- 
responding loop in the Grassmann manifold G(3;2), i.e., the 
space of two-dimensional subspaces of the three-dimensional 
state space spanned by |0), |1), and |e). By applying sequen- 
tially two TT pulse pairs with negligible temporal overlap, any 
desired holonomic one-qubit gate can be realized. 

We first consider the one-qubit case. A pair of zero- 
detuned laser pulses with the same real-valued pulse en- 
velope ^{t) couple selectively two ground state levels 
and 1 to an excited state e, as shown in Fig. [TJ For 
instance, by encoding 0, 1, and e in atomic or ionic hy- 
perfine levels, the laser pulses can be suitably polarized so 
as to assure that each transition is addressed separately. 
The corresponding Hamiltonian reads 

i/(i)(0 - r!(t)(L^o|e)(0|-fc^i|e)(l|+h.c.). (4) 

Here, the laser parameters cjq and wi satisfy |woP + 
|wip — 1, and describe the relative strength and rela- 
tive phase of the O e and 1 O e transitions. The 
Hamiltonian is turned on and off during the time inter- 
val [0,r], controlled by Vt{t). We take |0), |1) to define 
the one-qubit state space. 



A universal holonomic one-qubit gate can be realized in 
the above A system by choosing time-independent wq and 
LOi over the duration of the pulse pair. To see this, we first 
note that under this choice of laser parameters, the dark 
state \d) = — wijO) -I- wo|l) decouples from the dynam- 
ics, which in turn implies that the evolution is reduced 
to a simple Rabi oscillation between the bright state 
\b) = ujI\0)+ujI\1) and the excited state 0. The Rabi fre- 
quency is 0(i) and it immediately follows that the qubit 
subspace M(i) spanned by \ipk{t)) = e'^./'o ^"' (*')''*>) = 
U{t,Q)\k), k = 0,1, undergoes cychc evolution if the 
pulse-pair satisfies /„ Vt{t')dt' = tt. Secondly, the evo- 
lution is purely geometric since {'\lik[i)\H^^\t)\tpi[t)) = 
(A:|F(i)(i)|0 = for t e [0,t]. Under the above condi- 
tions, the final time evolution operator U{t, 0) projected 
onto the computational space spanned by {|0), |1)} de- 
fines the holonomic one-qubit gate 



C/W(C„) = n-cr. 



(5) 



where n is a unit vector and cr = {(Jx,cry,az) are the 
standard Pauli operators acting on |0),|1). By let- 
ting wq — sm{9/2)e^'^ and uji = —cos{9/2), we find 
n = (sin 6* cos 0, sin 6* sin 0, cos 6*). U^'^^Cn) is a univer- 
sal one-qubit gate. This can be seen explicitly by noting 
that two pairs of laser pulses corresponding to the unit 
vectors n and m applied sequentially results in 



= m ■ n + icr ■ (m x n). 



(6) 



This is an SU(2) transformation corresponding to a rota- 
tion of the qubit with an angle 2 arccos (n • m) around 
the normal of the plane in R'^ spanned by n and m. 
Here, Cn and Cm are loops in G{3; 2) based at M(0) 
and C = Cm o Cn- By suitable choices of n and m, any 
desired one-qubit gate can be realized. For instance, the 
choice n = (cos 0, sin (/), 0) and m = (cos0',sin0', 0) re- 
sults in the phase shift gate \k) ^ e^'''('^' '''') \k),k^ 0, 1, 
up to an unimportant overall phase factor. A Hadamard 
gate |fc) ^ -i=[(-l)'=|fc) + \k® 1)], fc = 0, 1, can be im- 
plemented by a single pulse with n = 4=(1, 0, 1). 

To describe the geometric nature of the above holon- 
omy, we lift the loop Cn to a loop Cn in 5(3; 2). As 
noted above, each such lift corresponds to a choice of 
gauge. The loop Cn in G{3; 2) may be represented by a 
set of complex 2-planes, each of which spanned by the 
single- valued vectors 



ICiW) = U{t,0)\d)^\d), 
IC2W) = e^^(*)C/(i,0)|5) 

= e'^'-''^ [cosS{t)\b) - ismS{t)\e)] 



(7) 



in the three-dimensional complex vector space C"'. Here, 
S(t) = J^n{t')dt' and the global phase factor e*''^*) has 
been inserted to ensure that ^2(1')) = |C2(0)). The loop 



Cn can be visualized by noting that |Ci) points in a fixed 
direction in C'^ around which IC2) rotates. Physically, |Ci) 
represents the dark state and IC2) describes Rabi oscilla- 
tions between the bright and excited states [9| . The oscil- 
lations correspond to a loop €„ in 5(3; 2) represented by 
the single-valued gauge choice in Eq. ([7]) that projects 
onto the loop Cn of complex 2-planes in G{3;2). The 
connection one-form associated with this gauge reads 




n{t)dt 

which results in the holonomy matrix 

1 
-1 



[/(C„) = Z = 



(8) 



(9) 



for a single Rabi oscillation. An explicit calculation con- 
firms that Ek,i Zki\Ck{0)){Ci{0)\ = n ■ a ^ C/W(C„). 
Similarly, for the composite path C = CmoCn, we obtain 

Efe,i[/H(C^)ia(0))(0(0)|=m.n+z<T.(mxn) = f/(i)(C) 
from the holonomy matrix 



U{C) = W'ZWZ. 



(10) 



Here, the unitary overlap matrix with components 
Wki = (Cfc(O)IO(O)) corresponds to an integration of a 
pure gauge connection one-form iV'dV along any path D 
in 5(3; 2) that connects the initial bases {|Ci(0))j I C2( 0))I 
and (ICUO)), |C2(0))} of C„ and C^, respectively [lfl|. In 
other words, the loop C = Cm ° Cn in G(3; 2) is rep- 
resented by the loop C ~ V^^ o Cm oV o Cn in 5(3; 2), 
where the four path segments correspond to the four non- 
commuting factors in the right-hand side of Eq. (jlOl) . 

To complete the universal set, we propose a physical re- 
alization of a non-adiabatic holonomic two-qubit gate in 
an ion trap setup. Our scheme is a non-adiabatic version 
of Ref. j2jj which utilizes the S0rensen-M0lmer setting 
111 ] to design a holonomic two-qubit gate. The system 
consists of an array of trapped ions, each of which ex- 
hibiting an internal three- level structure 0, 1, and e. The 
transitions <^> e and 1 ^ e for an ion pair in the ar- 
ray is addressed by lasers with dctunings ±1/ ± S and 
±i/ =p (5, respectively, where 1/ is a phonon frequency and 
S is an additional detuning. Off-resonant couplings to 
the singly excited states |eO), |0e), |el), and |le) can be 
suppressed by choosing the Rabi frequencies |r2o(OI ^-^id 
|rii(i)| smaller than u [ll|. In this way, the effective two- 
ion Hamiltonian in the Lamb-Dicke regime reads 



ij(2) = '-L{\no{t)\'ao{(t>)^aoW 

-|17i(t)|2ai(-0)®ai(-(/))) 



(11) 



Here, rj is the Lamb-Dicke parameter (77^ ^ 1) and 
crfc((/)) = e''^/\){k\+h.c., fc = 0, 1. Note that due to the 
non-adiabatic nature of our gate, the ancilla state a of the 
original adiabatic scheme in Ref. [2] is no longer needed. 



The phase and ratio |f7o(t)P/|f)i(i)P = tan(0/2) 
should be kept constant during each pulse pair. These 
requirements imply that the two-qubit Hamiltonian ma- 
trix H\,i^^{t) vanishes, which ensures that the evolu- 
tion is purely geometric. Under the tt pulse criterion 
^ J^ ^\na{t)\'^ + \ni{t)\'idt = TT, the final time evolu- 
tion operator U{t, 0) projected onto the computational 
space 1 00), |01), |10), |11) defines the two-qubit gate 



C/(2)(C) = cos6l|00)(00| + |01)(01| + |10)(10| 
-cos6l|ll)(ll|+sin6le"^'^|00)(ll| 
+ sin6ie*'*|ll)(00|. 



(12) 



C is a path in Q{3; 2) as the effective dynamics takes place 
in the three dimensional subspace spanned by |00), |11), 
and |ee) of the internal degrees of freedom of the ions. 
Due to its entangling nature, U^-^HC) is universal when 
assisted by one-qubit gates |12| . 

In practical implementations, the gates need to be ro- 
bust to the error caused by the finite life-time of the ex- 
cited state e. To test this, we add decay of e to the non- 
adiabatic holonomic one-qubit gate. We compare the re- 
sulting fidelity with that of the corresponding adiabatic 
gate. We consider the phase shift gate \k) — > e**'''^/2|fc), 
k — 0,1, in the non-adiabatic and adiabatic scenarios. 
This gate can be implemented adiabatically utilizing the 
A- type system, but now by varying the two laser cou- 
plings independently so as to remain approximately in an 
instantaneous dark state in the limit of large run-time T. 
We assume that the excited state decays to the auxiliary 
ground state level \g) with rate 7. We model the decay 
with the Lindblad equation 



et 



= -,-[ff(i) 



[H^'Ht), Qt] + 2LgtL^ ~ L^Lgt ~ QtL^L, (13) 



where gt is the density operator and L — y/^\g){e\. Fur- 
thermore, ff(i)(i) = n(t)(u;o\e){0\+uji\e){l\ -l-h.c.) and 
i?(i)(i) = n {uJi{t/T)\e){l\ + uJait/T)\e){a\ + h.c.) is the 
Hamiltonian in the non-adiabatic and adiabatic settings, 
respectively. In the adiabatic case, note that the 0- 
state is decoupled from the excited state and that the 
a state is another ancillary ground state level [2|. Two 
hyperbolic secant tt pulse pairs are chosen to imple- 
ment the non-adiabatic phase shift gate. Explicitly, we 
choose n{t){u}o,uji) = /3sech(/3t)(-l, 1)7^2 and n{t - 
At) iuj'o,uj[) = Psech[P{t - Ai)](-1, e'''''/'^)/V2, where /3 
is the amplitude of the pulses and At is the temporal sep- 
aration of the two pulse pairs. The ideal adiabatic gate is 
generated in the T — >■ cx) limit by varying the laser cou- 
plings oji = sin(^/2)e"^ and oja = — cos('i?/2) along the 
loop {d,ip) = (0,0) ^ (f ,0) ^ (f,7r) ^ (0,7r) ^ (0,0) 
at constant speed. 

In Fig. U we show the fidelity {^\W (C) goutU {C)\0 , 
computed numerically for 4000 input states |^), uni- 
formly distributed over the Bloch sphere. Here, U{C) is 
the non-adiabatic or adiabatic holonomic gate and gout is 



the output state computed from Eq. p3)) . The fidehties 
are plotted as functions of the dimensionless quantities 
/3/7 and J7T in the non-adiabatic and adiabatic cases, 
respectively. Note that the pulse duration in the non- 
adiabatic setting decreases with increasing /3/7 since the 
pulse area is set to the fixed value it. Thus, by increasing 



J 



P/j we effectively speed up the gate. Furthermore, we 
have chosen fl/j = 12.5 and jAt = 8, where the lat- 
ter choice guarantees that the pulse overlap is negligible 
for the /3/7 range shown in the figure; a necessary con- 
dition to avoid any spurious dynamical contributions to 
the gate. 




FIG. 2: Influence of decay with rate 7 of the excited state e on the non-adiabatic and adiabatic holonomic phase shift gate 



|fc) 



^ikw/2 



\k), fc = 0, 1. The effect is quantified in terms of minimum (blue), average (black), and maximum (red) fidelities. 



The three panels show from left to right the non-adiabatic gate with decay, the adiabatic gate with decay, and the adiabatic 
gate without decay. Choosing hyperbolic secant vr pulses with amplitude /?, the non-adiabatic fidelities are plotted as functions 
of the dimensionless quantity P/j. We plot the adiabatic fidelities as functions of the dimensionless quantity fiT, where Q 
is time-independent global strength of laser couplings and T is the run-time of the gate. We have chosen n/7 = 12.5 and 
"fAt = 8, where At is the temporal separation of the two laser pulses in the non-adiabatic setting. At is chosen sufficiently 
large to guarantee negligible pulse overlap for the /3/-y range shown in the left panel. 



r 



The fidelities of the non-adiabatic gate tend monoton- 
ically to unity in the large /3/7 limit (left panel). This 
demonstrates that the non-adiabatic version of the holo- 
nomic phase shift gate can be made robust to decay of 
the excited state by employing sufficiently short pulses. 
A key point with adiabatic holonomic quantum compu- 
tation is that the population of the decaying excited state 
becomes negligible in the adiabatic limit. This behavior 
is confirmed as the fidelities of the adiabatic gate in the 
presence of decay tend to unity in the large T limit (mid- 
dle panel). The oscillatory behavior, on the other hand, 
is due to non-adiabatic effects originating from the finite 
run-time of the gate and is thus present also when the 
decay is set to zero (right panel). The revivals seen in the 
fidelities of the adiabatic gate without decay have been 
pointed out previously in Ref. J13l |. 

In conclusion, we have developed a non-adiabatic gen- 
eralization of holonomic quantum computation (HQC) 
with the primary purpose to find ways to construct uni- 
versal sets of robust high-speed geometric quantum gates. 
We have demonstrated an explicit realization of a uni- 
versal set of holonomic one- and two-qubit gates in non- 
adiabatic evolution in three-level A configurations. The 
scheme requires coherent control of fewer levels and be- 
haves simpler under decay of the excited state compared 
to the holonomic gates proposed for adiabatic evolution 
in tripod configurations [2|-|4| . Our gate opens up for the 



possibility to perform universal quantum computation on 
short-lived qubits by purely geometric means. 
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